Switching the decay rate of an emitter inside a cavity in the time domain 
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We have theoretically studied the deterministic switching in the time domain of the spontaneous 
emission of two-level emitters in an optical microcavity. We consider a system in the weak coupling 
regime, for which the refractive index of the cavity material is modified on a time scale shorter 
than the emitted lifetime, through the injection of free charge carriers. The equation of motion 
is introduced for the excited population of two-level emitters coupled to a time-dependent local 
density of optical states. We then derived the emitted intensity as a function of time for both 
continuous wave and pulsed excitation of the embedded emitters. We demonstrate that a short 
temporal increase of the radiative decay rate drastically increases the emission intensity during the 
switch time for both continuous wave excitation and pulsed excitation. For pulsed excitation, the 
resulting spontaneous emission shows a distribution of photon arrival times with a deterministic 
burst of photons during the switch event, that strongly deviates from the standard exponential law. 
Finally, we derive a figure of merit that quantifies the total switching action. 
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I. INTRODUCTION 

Light is essential for myriad processes around us: in nature, to human life, and to technological applications and 
everyday appliances. It is known that an elementary light source such as an excited atom, quantum dot or molecule 
emits a photon either spontaneously or stimulated by an external field [l[ . The rate of spontaneous emission is not 
an immutable property of a light source [H, 0| but the rate strongly depends on its surroundings. For an emitter 
coupled to a continuum of electromagnetic modes, quantum mechanics predicts the rate of spontaneous emission of 
an excited two-level source to follow Fermi's golden rule [|[: the rate is determined by a product of atomic transition 
matrix elements of the dipole operator with the local radiative density of optical states (LDOS), that typifies the 
surroundings [f|. The LDOS counts the number of modes in which a photon can be emitted, and it can be interpreted 
as the density of vacuum fluctuations at the position of the source. A main feature of spontaneous emission is its 
dynamics: an emitted photon is detected at a random moment in time after the emitter has been excited. Both the 
distributions of emitted photons and of the excited-state population decay exponentially in time. 

A well-known tool to modify the average spontaneous decay rate of a source in the frequency domain is a resonant 
cavity tuned to the source's emission frequency. It was first realized by Purccll that the decay rate of an atom can be 
increased @, known as the Purcell effect. In the weak coupling regime the Purcell factor F p gauges the change of the 
radiative decay rate r cav ; ty of an emitter in a cavity relative to the rate To in a homogeneous medium: F p = r cav ;ty/ro. 
Following the pioneering work by Gerard et al. Ml, many groups have demonstrated the Purcell effect with quantum 
dots embedded in solid-state microcavities [8j-Fl0j| . To date, impressive progress has been achieved in controlling 
spontaneous emission in the frequency domain with nanophotonic structures [ll| - [l3| , like microcavities and photonic 
crystals 0, [H and nano-antennas [161 ]. In all cases, however, the modification of the LDOS is stationary in time. 
Thus, the radiative decay rate is time independent and the distribution of photon emission times is single exponential 
in time and completely determined by this rate. 

In this work, we theoretically propose a novel paradigm in light emission by modifying the environment of an 
ensemble of two-level sources in time during their lifetime. By utilizing fast optical switching of a microcavity, we can 
tune the cavity resonance and drastically change the LDOS at the emission frequency within the emission lifetime fvj\ . 
As a result, we anticipate bursts of dramatically enhanced emission, concentrated within short time intervals. The 
spontaneous emission process remains stochastic but results in a strongly non-exponential temporal distribution that 
is comple tely controlled by the experimentalist. Our approach offers a tool to dynamically control the light-matter 
coupling 1 1 81 ] and, in the long run, offers a novel tool to realize the breakdown of the weak-coupling approximation 
and realize non-Markovian dynamics in cavity quantum electrodynamics. 

First we introduce the equation of motion that describes the population density of an ensemble of two-level sources 
with a time-dependent decay rate due to a time-dependent LDOS. From this equation, we derive the time-dependent 
emitted intensity for continuous wave and pulsed excitation of an ensemble of two-level sources such as quantum dots 
[3] or rare earth atoms [l9j in a cavity. Since micropillar cavities are known to be a versatile class of microcavities we 
choose them as an example. The decay rate of the ensemble, determined by the LDOS, is switched by exciting free 
carriers, which is a well-known control mechanism in the time domain for nano-cavities [20l - [25| . 

II. SWITCHING THE DECAY RATE 
A. Rate equations 

We consider a two-level emitter in a medium with a strongly dispersive LDOS p(w,r) in a photonic microcavity 
and investigate what happens if the LDOS changes in time, thereby modifying the radiative decay rate in time. To 
derive the rate equation of a two-level source we start with the equation of motion of the probability amplitude of the 
excited two-level source c a (t) [2(| with a LDOS p(u),ed,r,t') that depends on time t' 

^dT = -^/7o 00cQ( '' V/9(w ' ed ' r ' 0e4(aJ " aJd)(t '^ )M '- (1) 

Here d and are the amplitude and orientation vector of the transition dipole moment, respectively, h the reduced 
Planck's constant, eo the dielectric constant of vacuum, r the emitter position, and u>d the emission frequency. 
For convenience, we only write the time dependency of c a (t), but it should be kept in mind that the amplitude 
c a {t, r, e c i,LUd) also depends on r, and iOd [271 ] . 

In the following we limit ourselves to the weak coupling limit where the factor (wp(w, e<j, r, t')) varies slowly over 
the linewidth of the emitter. This approximation is known as the Markov approximation [28j or the Wigner- Weisskopf 
approximation [l[. In the Markov approximation we can take uip(w, e^, r, t') out of the frequency integral and Eq. JT]) 
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can be simplified to [29( 

dc a (t) d 2 



dt 2heo Jo 

The integral in Eq. © can be evaluated to yield [3(| 

dc a (t) d 2 



Ca(f)ir5(t - t')uj d p(u} d , e d , r, t')dt'. (2) 



which can be written as 



with r rac j the radiative rate 



dt 2heo 

dc a {t) _ r rad 
dt ~ 2 



c a (t)TTUJdp(u}d,ed,r,t), (3) 



Ca(t), (4) 



d 2 u d 'K 

r ra d(i) = — r p[uJd,e d ,r,t). (5) 

tieo 

Equation ([5|) is Fermi's golden rule [H augmented with a time-depend LDOS. This shows that in the Markov limit the 
instantaneous radiative rate r rad (i) directly follows the time dependence of the LDOS. In case of a time-independent 
LDOS the rate r rad (t) = r rad is constant and Eq. §4§ shows the well-known feature that the amplitude c a (t) decreases 
exponentially with the rate Similarly, the probability |c a (i)| 2 of the two- level system to be excited decreases 

exponentially according to 

\ Ca (t)\ 2 = MO)! 2 ^- 4 . (6) 

For a time-dependent LDOS the rate in Eq. ([6]) is no longer constant and the excited state population decreases 
non-exponentially and thus deviates from the standard Markovian dynamics. 

With Eq. ((4]) we can write the equation of motion for the population density N 2 (t) for an ensemble of N identical 
(no inhomogencous broadening) and non-interacting two-level sources. To complete the model we include a time- 
dependent excitation term for the sources and a non-radiative decay rate r nra( j. The equation of motion for the 
population density becomes 

dN 2 {t) P mc (t) 

77 = ?7abs-£ (r ra d(i) + r, lrad ) N 2 (t). (7) 

dt flLO e xc 

The first term describes excitation and depends on the excitation power P exc {t), the excitation frequency uj exc and, 
the absorption efficiency of the excitation power that reaches the two- level source r^abs ■ The second describes radiative 
decay and the third term non-radiative decay. For convenience, we write iV 2 (t) only as function of time t in Eq. ([7J . 
For any sub-ensemble of the N two-level sources A^i) also depends on r, ed and oj d , hence for an inhomogeneously 
broadened (sub-)ensemble we should average over r and e d - At any rate the general solution of Eq. ([7]) is 

N 2 (t) = N 2 (Q) + f (ij^^^l - (r rad (f) + r nrad ) N 2 (t')) dt'. (8) 



The corresponding radiated emission intensity I em (t) is given by [31( 

Iem(t) =T Iad (t)N 2 (t), (9) 

which means that the total emitted light intensity is proportional to the instantaneous radiative decay rate and the 
population density. Equations (jSJ and © form the basis for our further discussion and they will be used to calculate 
the emission for an ensemble of emitters with a time-dependent radiative rate in the cases of continuous-wave and 
pulsed excitation. 



B. Time dependent radiative decay rate 



The central goal of this work is to describe the effects of a time-dependent radiative decay rate r rad (t) that is 
realized by dynamically changing the LDOS in time at the position and frequency of an emitter. In general we 
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FIG. 1. (color online) Radiative decay rate normalized to the unswitched rate To (solid line) as a function of time after exciting 
the emitter. The two thick curves show the result of a switching event at to pu — 10 ps that either enhances (long dashed) or 
inhibits (short dashed) the decay rate by a factor of 5. The modified decay rate relaxes back to the unswitched rate within the 
effective switching time of r sw = 35 ps after the switching event. 



can separate the time-dependent decay rate into a constant rate Tq and a time-dependent change in the decay rate 

Ar rad (t) 



rad 

(t) = r + Ar rad (*). 



where A ra d(i) is proportional to the change in the LDOS 



Ar rad * = ——*Ap(t). 

he 



(10) 



(11) 



We assume that the time-depended part is the result of a short switching event that quickly changes the LDOS within 
a characteristic switching time t sw . Two examples of normalized time-dependent decay rates are shown in Fig. [1] for 
emitters excited at t = ps followed by a short switching event at t = t 0pu = 10 ps. In the examples in Fig. Q] the 
switching gives either an enhancement or inhibition by a factor of 5 that decay with a characteristic effective switching 
time of r sw . 

In the following we choose for practical reasons a scheme where the emitter is embedded in a semiconductor cavity. 
Moreover, the LDOS is modified in time by controlling the refractive index by means of the free carrier density in the 
semiconductor, as excited by a short optical (or electrical) pump pulse at t = t pu . The induced change in the refractive 
index is proportional to the free carrier density (32[ and the resulting change in the LDOS depends strongly on the kind 
of the microcavity [33j . The excited free carriers recombine exponentially with a characteristic recombination time 
r sw , after which the refractive index is restored to its original value [HI [HI]. Here we use r sw = 35 ps, characteristics 
for GaAs [p. 

As an illustration, we consider in Fig.[5]the effect of switching the resonance frequency w caVj o of a microcavity with a 
Lorentzian LDOS with linewidth r cav in the spectral vicinity of the emission frequency u>d of a emitter. The decrease 
in the refractive index induced by the free carriers will lead to a positive frequency shift of the resonance frequency 
Wcav(i) as indicated in Fig. [5J The emitter is initially detuned from the cavity resonance and experiences a low 
radiative rate Torad. During the switch event the cavity peak is tuned into resonance with the emitter as shown as the 
dashed Lorentzian in Fig. [5] This change results in a rapid increase in the LDOS at the emitter frequency and greatly 
enhances the decay rate r ra d(i) from the initial value r ra d(0) = Tq to its maximum value of r ra( j(A<) = Tq + Ar rac i 
and back to Tq at a time At. The effective switching time in this scenario is therefore given by 



At 



l^cav 

(At) — w cav ,o| 



r« 



(12) 



A shorter effective switching time can thus be realized by either faster tuning of the cavity resonance in a time A or 
by increasing the spectral tuning range relative to the cavity linewidth within the time At. 

We note that this switching procedure is very flexible and we can effectively move along different trajectories on the 
cavity's LDOS by choosing the initial detuning and strength of the switching effect. For example, we can consider an 
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FIG. 2. (color online) Schematic graph of the switching process as experienced by a quantum emitter (green) emitting at a 
frequency Ud in the spectral vicinity of a cavity resonance whose frequency is switched in time. The cavity has a Lorentzian local 
density of states (solid line). Initially the emitter is detuned from the cavity o by nearly one cavity linewidth, 

leading to an effective radiative rate To. The switching process moves the cavity resonance up in frequency w ca v(£) (gray dashed). 
The cavity is then tuned into resonance with the emitter that thus experiences a decay rate strongly enhanced by Ar ra( i. Within 
one cavity linewidth from the resonance, switching of the cavity resonance can be approximated by a linear shift of the decay 
rate versus frequency (red dashed line). 



alternative experiment where the emitter starts on resonance with a radiative rate that is already Purcell enhanced. 
The switch then detunes the cavity resonance away from the emitter's frequency and thus inhibits the spontaneous 
decay rate. The steep slope of the cavity LDOS gives a rapid change in the LDOS that can be used to either greatly 
enhance or inhibit the radiative decay rate, relative to the unswitched rate. 

For initial dctunings between the cavity and emitter frequency smaller than the cavity linewidth (wd — ^cav.o < r cav ) 
we can approximate the steep slope of the Lorentzian resonance as a linear trend shown as the red dashed line. Wc 
can therefore effectively approximate a linear relation between the excited free carrier density and the radiative decay 
rate. For a typical switching pulse with a Gaussian temporal width r pu = 120 fs, much shorter than the carrier 
recombination time (35 ps, see Ref. [2214241 ] ) . we can separate the excitation and relaxation time scales of the free 
carriers. Using the linear relation between the carrier density and the decay rate discussed above, the time-dependent 
decay rate 



r ra d(£) = r + Ar rad e _ 



Q(t - t 0pu ,T pu ) 



(13) 



is then decomposed into an constant decay rate r and a change induced by the switch that is turned on at time 
^Opu- The magnitude of the switched term in Eq. (|13[) then decays exponentially with an effective switching time 
comparable to the free carrier relaxation time. The initialization of the switch is evaluated as an integral over the 
Gaussian switch pulse resulting in an error function 



0(t,r) 



2 1 

7T T 



/ e 2t- 



dt' = - 



1 + erf 



V2t 



(14) 



that in the limit r pu -C t sw can be approximated with a Heaviside step function. 

The normalized time-dependent decay rate in Fig. [T] giving an enhancement thus depicts the situation where we 
quickly tune the cavity frequency, initially off-resonance, into resonance with the emitter as illustrated in Fig. [2J As 
a result the radiative rate is quickly and greatly increased at t = 10 ps before decreasing at a speed set by switching 
time. Similarly, the other curve illustrates the situation where the light source is initially on resonance and the cavity 
is then switched out of resonance. In the examples in Fig.[T]thc switching gives either an enhancement or inhibition by 
a factor of 5, which is a realistic change observed on ensemble of quantum dots in micropillar cavities Q- Note that a 
constant relative decay rate of r ra d(i)/Fo = 1 corresponds to the unswitched case, typical for all Purcell experiments 
performed to date pTHlol [33l| . Most striking is the fast dynamics in the decay rate: both the switch pulse duration 
T pu and the exponential decrease with decay time t sw are much faster than the intrinsic lifetime I/Tq = 1 ns typical 
for quantum dot emitters. 

With the time-dependent decay rate given in Eq. (fT3")) we can state that the assumptions for the population 
dynamics (from Eq. (Q]) to Eq. ([!]))) arc valid if uJdT pu 1. This is the case for a pump-pulse duration r pu = 120 fs 
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FIG. 3. (color online) (a) Population density and (b) emitted intensity normalized to the steady state values for a CW-excited 
ensemble of emitters, such as quantum dots, for a free carrier switch event at t — t pu = ps. The green dashed curve shows the 
influence of a change in the decay rate for a Purcell factor of 5 (Ar ra d = 4Fo and To = 1 ns~ ) and switching time identical to 
the free carrier relaxation time r sw = 35 ps. The red dashed line shows the residual quantum dot excitation by the free carriers 
with a power of P pu r] rcc = 0.5P Ocxc . 



and ujd = 2 x 10 s -1 (corresponding to = 950 nm). In this case the fastest variation of the LDOS in time has a 
characteristic time that is much longer than the width of the delta function resulting from the exponential in Eq. ([1} . 

In the following sections we will focus on the effects of the time-dependent decay rate on the population density 
and emitted intensity of an ensemble of two-level sources. We will look at two examples: one where we switch an 
ensemble of sources that are continuously excited and one where the sources are excited by a short pulse (see Fig. [5]) 
and the environment is subsequently switched. 



C. Population and emission dynamics for continuous wave excitation 

In the case of continuous wave (CW) excitation of the two-level sources the direct excitation is constant in time 
Pexc(t) = Poexc = ^Wexc/^abs ■ We account for the possibility that a small fraction of the free-carriers used in the 
switching process will excite the emitters as discussed in Appendix [X] We therefore additionally include an excitation 
term that is proportional to the free carrier concentration. Substituting the time-dependent decay rate (Eq. (|13|l ) and 
excitation power (Eq. (|A1[) 1 into Eq. (|SJ) and integrating over time yields the population dynamics 

m) ^ + rya^exc t+ ^f f p Zi^d \ ^ 

nhJexc Jo V / 

-J ( r + Ar rad e ^ e(t - t pu,T pu ) + r nrad J N 2 (t")dt". (15) 

The first term Nq is the initial population density at time zero set to the steady state population density without 
the switch. The two next terms represent the pumping of the emitters and the last is the relaxation term from 
spontaneous emission. Since Eq. (|15|) cannot be solved analytically, we have solved it numerically and the results are 
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plotted in Fig. [3] for three examples. In the unswitched case with a continuous excitation and a time-independent 
radiative decay rate, both the emission intensity and excited state population are constant. This value is determined 
by the equilibrium between the spontaneous emission and excitation of the emitters. 

Figure (3}d shows the dramatic influence on the emission dynamics for a switch pulse at t = t pu = ps. The 
switching parameters are the same as in Fig. [T] with To = 1 ns, Ar rac j = 4Fo, ?7 rC c = and t sw = 35 ps. In the 
first 100 ps after the switch at t = 0, we see a strong burst of spontaneous emission with an intensity that is 6x 
higher than the steady state value. The switch therefore creates a short intense pulse above a weak background. The 
duration of the burst is much shorter than both the unswitched and switched radiative decay times of 1 /Tq = 1 ns 
and l/(r + Ar ra d) = 200 ps. respectively. The temporal width is defined by the switching time of 35 ps. After the 
switch event is completed at t ss 200 ps, the emission intensity decreases below the unswitched emission level. This 
drop is caused by the rapid depletion of the population density during the emission burst shown in Fig. [3^,. During 
this burst, the emission rate is indeed much larger than the (constant) pumping rate. After the switching event, 
the depleted population of excited emitters emits therefore less light than in the steady state regime. The recovery 
dynamics of the population of excited emitters and of the emitted signal follow an exponential law, governed by the 
radiative lifetime Tq. 

We next analyze the effect of a potential excitation of the emitters by the free carriers excited by the switching 
pulse. In Fig. [3] we show the population dynamics and emission intensity for r rac j = Tq and P pu ?/rcc = 0.5Pocxc- For 
these parameters the switching pulse does not change the decay rate but only excites the emitters. We therefore see 
a small increase in the population density just after the switching pulse. This increase then decays back to the steady 
state level with a time constant given by the radiative decay time 1/Fo. The corresponding emission dynamics shows 
an almost negligible increase in intensity. Even when the free carrier excitation rate equals the optical free carrier 
excitation rate (-Pp U 7?rec = -Fbexc) the emission intensity is only enhanced by 3%. The small increase is caused by the 
short duration of the pump. Therefore the difference in time scales between the carrier relaxation time r sw and the 
radiative decay time 1 /r ra d (t) governing the intensity dynamics let us distinguish between true decay rate switching 
or emission increase due to residual free carrier excitation. 

From Fig. [3] we conclude that switching the decay rate leads to a strong enhancement of the emission intensity 
compared to the unswitched case. The duration of the emission enhancement is determined by effective switching 
time, while the timing of the emission burst is determined by the arrival time of the switching pulse. By applying 
the switch pulse we are able to detcrministically enhance an otherwise small and constant radiative decay rate and 
to create an emission burst, whose duration is only limited by the effective switching time. 



D. Figure of merit for pulsed excitation 



In this section we study the dynamics for the excited population of emitters after a pulsed excitation, when the 
environment is subsequently switched during their decay time. It will appear that the dynamics for pulsed excitation 
strongly differs from the dynamics under continuous wave excitation. Wc assume that a short excitation pulse with 
amplitude Poexc initializes the system at t — <ocxc such that we have an initial population density N%(t = toexe) = N02. 
At times after the excitation pulse the dynamics of the population density is governed only by the time-dependent 
decay rate and this gives a monotonous decrease in the population density. If we approximate the short excitation 
pulse by a Dirac delta pulse P CX c(t) = S(t — iooxc)^0cxc in the rate equation (Eq. (0) it can be solved analytically for 
times after the excitation (t > iocxc)- In this case Eq. (|7|) simplifies to 

— — ° eXC ^ = — (r rac i(i — iooxc) + r nra d) N 2 (t — tooxc), (16) 

at 

which can be integrated to yield 

N 2 (t-t Qexc )=N (a exp(j -(r rad (i')+r nrad )dt'^ . (17) 

Equation 1171 describes the population density for any time-depended decay rate F ra d(i) as a function of time t after 
the excitation process is over. Inserting the switched decay rate Eq. (| X3[) into Eq. p7|) and solving the integral over 
the constant part of the decay rate yields 

N 2 (t - Wc) = Ar 02e -( r »+ r °"^)(*-*°~c)-Aa rad (t) (18) 
where we have defined a dimensionlcss time-dependent switch parameter Aa ra d(£) 

Acw(t)= / Ar rad e (t ^w Pu) Q(t - t 0pu , r pu )dt'. (19) 
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This parameter is a figure of merit that describes the change in population density due to the change in the decay 
rate. A negative Aa rad (t) results in a population density that decays slower compared to the unswitched situation, 
while a positive Aa ra{ j(t) results in a faster decay. If we assume that the duration of the switch pulse r pu is short 
compared to the effective switch time r sw , the integral in Eq. (|19p can be split into two parts - before and after the 
switch t = r pu - and Aa ra d simplifies to 

/ -U-*0pu) \ 

Aa rad (i) = Ar rad T sw I 1 - e -™ ) 6(t - t 0pu , r pu ). (20) 

Here Q(t — topu, T pu ) is a step function from to 1 that accounts for the fact that there is no change in the decay rate 
before the switching pulse arrives at t = iopu- In the limit of time t going to infinity Aa rad (i) becomes 

Aaoo = lim Aa rad (i) = Ar rad r sw . (21) 

t— >oo 

Equation (|2ip shows that Aa ra( j(t) is nonzero even in the long time limit and is given by a product of the magnitude 
Ar rad of the switch and the effective switch duration. The switch therefore has an effect on the population dynamics 
even long after the switch event. The dimensionless switching parameter Aaoo is therefore a useful figure of merit for 
the total switching effect on the excited state population. 

We can quantify the effect the switch has on the population at long times by using Eq. (|2"Tj) and the limit in Eg. (fT5)) 
to calculate the ratio between the switched population density N 2s and the unswitched population density A^us m 
the limit of t going to infinity 

lim J^> = e -Arw 3W = (22) 
t^°°N 2us (t) y ' 

Equation (|22j) quantifies the long term effect of the switching on the population density as a result of a momentarily 
short change in the decay rate. 



E. Population dynamics for pulsed excitation 

Figure |4ji displays the excited state population Eq. (|18|) for three cases: Without switch pulse the decay rate is 
stationary and as expected the population and intensity decay exponentially with the initial rates Tq. Secondly the 
green long dashed curve shows a switch that tunes a cavity into resonance with an emitter and induces an enhanced 
decay rate by a factor of 5 (Ar rad = 4To) from To = 1 ns~ J . Thirdly the red short dashed curve represents the 
opposite case where a cavity is tuned out of resonance by the switch and induces an inhibition in the decay rate 
by a factor of 5 starting from a high initial rate To = 5 ns -1 . For the two switched examples (dashed lines) the 
population density clearly decays non-exponentially. In the enhanced case the population decays exponentially before 
the switching pulse with the same rate as in the unswitched case. As the switching pulse arrives at t = t pu = 150 ps 
the population decreases faster and thus deviates from exponential decay. During the effective switching time of 35 ps 
the population density continues to deviate from an exponential decay. A few recombination times after the switch 
event (t > 250 ps) the decay returns to its original rate but the absolute value of the populations is reduced compared 
to the unswitched case. Using Eq. (|T8| and the figure of merit (Eq. pip ) we see that the larger decay rate induced 
by the switch depletes the excited state population faster, thereby lowering the population density at long times. 
The population density does not converge back to the unswitched curve (unlike the continuous wave case) since there 
is no continuous reexcitation of the emitters. The situation is reversed for a switch that induces an inhibition of 
the spontaneous emission : the population experiences also a non-exponential decay after the switch; however, the 
population is now larger than its reference value (unsswitched case) at long times. 



F. Emission dynamics for pulsed excitation 



In this section we discuss the corresponding effects on the dynamics of the emitted intensity from emitters in a 
switched environment. According to Eq. ([9]) the emitted intensity I{t) is given by the product between the excited 
state population and the radiative rate that is also time-dependent. Modifications to the decay rate are therefore 
directly reflected in the total intensity emitted by the emitter. So, for a large dynamic change in the decay rate, we 
also expect large changes in the emitted intensity. One striking consequence is that for a time-dependent decay rate 
the population density and the emission intensity are no longer directly proportional, contrary to the steady state 
case 1311. 
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FIG. 4. (color online) (a) Time resolved population density for an emitter excited at t — t exc — ps showing the effect of 
two different switch events at t = t pu — 150 ps. The chosen parameters model the effect of a switch event that either tunes a 
cavity resonance into (green long dashed) or out of (red short dashed) resonance with the emitter frequency. Without switch 
the populations decay exponentially with a rate of To = 1 ns _1 and To = 5 ns — , respectively, in the two examples (solid lines). 
The switch event leads to an enhanced or inhibited radiative rate by a factor of 5 relative to the unswitched rate. These time- 
dependent rates result in a short decrease or increase in the populations relative to the unswitched cases. At long times after the 
effective switching time r sw = 35 ps, the slopes tend to their initial values for both examples, (b) The corresponding spontaneous 
emission intensities from the quantum dot relative to the initial values after excitation for the same two examples presented in 
(a). The small changes in the population density manifest themselves as large changes in the emitted intensity. Switching the 
cavity into resonance with the quantum dot (green long dashed) results in a sharp burst of intensity with a temporal duration of 
Tsw Tuning out of resonance leads to a fast drop in the intensity. 



Inserting the dynamic decay rate Eq. (flU)) and the population density Eq. (fl~5|) into Eq. © with r) tec = yields 

i(t) = (r + Ar rad (t)) Ar O2e - r o tot (t-toe«)-A« rad ( t)0 ^ _ fQexc) Toxc)j (23) 

where Ar ra d(<) = Ar rac ie t ™ 0(f — t 0pu , r pu ) , and Aa(t) is given by Eq. (f^Uj) . The main difference between the 
population density dynamics and the emitted intensity is the presence of the decay rate prefactor (ro + Ar ra d(i)). 
Nevertheless, the intensity in Eq. (J23J) is still proportional to the population density so that the influence of the 
switching process will remain visible in the emission intensity even long after the switch event has passed as discussed 
in section Hi Dl The relative intensity to the unswitched intensity at long times is thus given by lim^oo I(t)/I us (t) = 
e -Ar rad T sw ag ^ ne ex p 0nen t is the same as in Eq. §TZ\> and the time-dependent decay rate Ar rac i(t) in the prefactor 
tends to zero for long times. 

Figure |4Jd shows the normalized time-resolved emission dynamics for the same two cases as for the population 
density in Fig. 0^: one where the radiative rate is quickly enhanced from an initial low rate of Tq = 1 ns _1 and 
another where the radiative rate is inhibited from a high value of Fo = 5 ns _1 . The emitters are first excited at 
t ex = ps and the emission intensities then start to decay exponentially with the same rate as the population density 
as expected in the weak coupling limit. A switching pulse arrives at t = < pu = 150 ps whose effect is to either quickly 
enhance (green long dashed) or inhibit (red short dashed) the radiative decay rate from the initial rate by a factor of 
5. 

In the case where the switching pulse quickly enhances the decay rate we see in Fig. 2Jd a short and high intensity 



10 



burst on top of the normally decaying signal; the intensity thus strongly deviates from an exponential decay. The 
relative magnitude of the enhancement is equal to the maximum Purcell enhancement and the temporal shape closely 
follows the modulation in the decay rate as expected from Eq. ([23)) . In this example the temporal shape follows the 
exponential change in the decay rate and the temporal width is limited by the effective switching time of 35 ps. This 
time is much shorter than the minimum Purcell enhanced decay time of 200 ps. Let us note at this stage that the 
effective switching time could be further engineered to be as short as 2 — 3 ps by either decreasing the free carrier 
lifetime 3^ or increasing the frequency shift of the cavity resonance. For times much longer than the switch time, 
we see a lower intensity relative to the unswitched case due to the depletion of the population density of the emitter 
discussed in Sec. Ill Dt 

In the second case in Fig. HJd the switching event inhibits the decay rate, which results in a short temporal drop in 
the emission intensity relative to the stationary case as seen in Fig. HJs. The temporal duration of the drop is again 
limited by the switch time. The drop in intensity exemplifies a highly unusual shape for a decay curve where the 
radius of curvature, during the free switching time, is negative during the decay. In the traditional paradigm of steady 
state spontaneous emission such a negative radius of curvature would be unphysical. An exponential decay with a 
stationary decay rate or even a sum of exponentials with stationary rates will always result in a positive radius of 
curvature in the resulting decay curve. This shows the flexibility of the switching mechanism to shape the temporal 
emission distribution of the emitted photons at will. 



III. DISCUSSION 



Spontaneous emission can be described as a stochastic process for a photon to be emitted from an excited light 
source [l[ . It is thus impossible to predict the exact time when an excited source will emit a photon. On the other 
hand, the distribution of photon emission times, averaged over many excitation cycles, is usually well-known. In the 
weak coupling limit the distribution of photon arrival times is exponential and characterized by a single rate given by 
Einstein's A coefficient. Much interest has been devoted to controlling spontaneous emission by modifying this rate 
using the Purcell effect by embedding emitters in a nano-structured environment @, HI, fiol HH . Several schemes have 
been implemented to tune the decay rate in time such as gas deposition, temperature and, electronic gating. However, 
the modification in the rate has in all cases remained constant in time during a single decay process. For this reason, 
the distribution of photon arrival times thus remains single exponential, and is simply described by one enhanced 
or inhibited single exponent. We have here introduced a new paradigm with dynamic control of the spontaneous 
decay rate in time using all optical switching. The spontaneous emission process remains stochastic in time but the 
dynamical change in the decay rate results in a strongly non-exponential temporal distribution of photon emission 
times. The active switching process allows us to dcterministically control the photon distribution in time. We have 
shown (in Fig. [4]) that photon arrival times can be bunched in short bursts where timing and duration of the burst 
can be fully controlled by the experimentalist. Naturally, within these short emission pulses, the individual photons 
still arrive at unpredictable moments in time. 

On a more fundamental level spontaneous emission arises from the interaction between a single quantum emitter 
and fluctuations in the vacuum field at the emitter position. By dynamically modifying the environment of the emitter 
our approach gives direct temporal control of the local strength of the vacuum field on timescales much shorter than 
the excited state lifetime. As shown in Fig. @] this allows to manipulate the excited state probability for a quantum 
emitter in time and subsequently control the time dependence of the single photon wave function of the emitted 
photon. Such control opens great prospects in quantum information processing and allows to shape the photon wave 
function emitted by single photon sources, for example for optimal mode matching of photons |35j and to enhance 
the absorption of single photons [36], [13] ■ More generally, by dynamically tuning a cavity in the vicinity of a quantum 
emitter we can drastically modulate the light matter coupling between the emitter and the cavity mode. This offers 
interesting prospects where a system is modulated between the weak and strong coupling regime while emitting a 
single photon. 

For very fast switching events the decay rate of the emitter can no longer adiabatically follow changes in the 
environment and the decay rate is not proportional to the instantaneous LDOS but depends also upon the past 
history of changes in the LDOS. So even a coupled cavity-emitter system, that in the steady state case would 
be weakly coupled, can by fast switching of the environment be brought into a regime where the weak-coupling 
approximation breaks down. Our method thus offers a novel tool to realize non-Markovian dynamics in cavity 
quantum electrodynamics, namely by very fast modulation. 

For a large ensemble of emitters our approach offers a tool to implement a bright ultra-fast light source based on 
spontaneous emission with a low temporal coherence. This source has potentially much shorter pulse duration than 
electronically controlled LEDs. The photon statistics of such a source differs significantly from known laser action 
such as Q-switching or cavity dumping). An ultra fast low coherence source may find applications in speckle-free 
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FIG. 5. The time- dependent excitation power of an ensemble of quantum dots showing both the direct excitation at ioexc = 5 ps 
with a duration of 1 ps and the secondary excitation due to free carriers excited by the switch pulse at to pu = 10 ps. The free 
carrier density decays with the carrier recombination time o/r sw = 35 ps. 



imaging which requires low coherence [36 



IV. CONCLUSION 

We have studied the deterministic switching of the Purcell factor and thereby the decay rate for spontaneous 
emission of light sources in a cavity. We have introduced a model geared toward experimental validation using free 
carrier switching of micropillars cavities with embedded quantum dots. We have demonstrated that by dynamically 
controlling the radiative decay rate of emitters during the decay time the emission intensity can be drastically modified 
during the switch time for both continuous wave excitation and pulsed excitation. For pulsed excitation the dynamic 
decay rate reveals a strongly non-exponential distribution of photon arrival times. Finally, we have introduced a figure 
of merit of spontaneous emission switching, equal to the product between the radiative decay rate and the switch 
duration. 



Appendix A: Direct excitation and free carrier excitation 



In this appendix we elaborate on the excitation of the emitters, described by the term P exc {t) in Eq. ((5]). We 
consider that the carriers excited by the switching pulse may not be perfectly shielded from the quantum dot region. 
Thus, there are two main contributions to the excitation power. The first term results from the direct excitation 
pulse, which excites the ensemble of quantum dots sources. The time behaviour of the excitation pulse is modelled 
with a Gaussian function [39j . The second contribution results from the leakage of free carriers into the quantum dots 
created by the switch pulse. The contribution is therefore proportional to the free carrier density and has the same 
temporal behaviour as the free carriers. The total excitation power is thus given by 



Popuffre 



-C-'Opu) 



_ ^Opu, Tpu)- 



(Al) 



Here Pocxc is the amplitude of the excitation pulse, ioexc the time at which the excitation pulse excites the quantum 
dots, r oxc the duration of the excitation pulse and, P pu the amplitude of the pump pulse. The rest of the parameters 
are defined in relation to Eq. (fT3|) . 

An example of the time dependence excitation power from eq. (|A1[) is shown in Figure [5] first the short peak from 
the intentional excitation pulse at tooxc = 5 ps used to excited the quantum dots followed a longer peak resulting from 
leakage of free carriers generated by the switching pulse arriving at iopu = 10 ps. The magnitude of the excitation 
due to free carriers is taken to be -Fbpu^rcc = O.lPooxc- Since the leakage rate is proportional to the free concentration 
this contribution has a sharp increase at iopu and then decreases exponentially with the free carrier recombination 
time taken to be t sw = 35 ps. In a practical experiment it is desirable to reduce rj voc by shielding the quantum dots 
from the free carriers generated by the switching pulse. One way is the use of a cavity where the optical mode for 
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the switching cavity extends over two dielectric regions separated by an air gap. Switching the refractive index of 
one region of the cavity would change the combined cavity frequency whereas the quantum dots in the other region 
would be isolated from the excited free carriers. 
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